(1)
摘 要：咸海萎缩导致的生态危机是"一带一路"沿线面临的重大环境挑战，其中土壤风蚀是引发盐尘暴与土地退化的关键过程。为系统揭示中亚咸海及周边区域土壤风蚀的长期演变规律，本研究基于修订版风蚀方程（RWEQ），在Google Earth Engine（GEE）平台中融合多源地理空间数据，构建了1990–2020年中亚咸海及周边区域风蚀模量年际数据集。数据集整合了MODIS植被指数数据、ESA-CCI土地覆盖数据、SRTM高程资料以及ERA5-Land再分析气象数据，经统一坐标系统、分辨率重建与时空一致性处理后，生成空间分辨率为1 km、时间跨度为31年的风蚀模量序列。数据分析表明，咸海干涸湖床为区域风蚀的主要源区，且自2012年起风蚀核心区呈现由湖床内部向外围扩展的显著趋势，反映出侵蚀格局的动态演变特征。为验证数据集可靠性，将模拟结果与咸海沙尘监测站（2000–2005年）记录的粉尘通量实测数据、气溶胶光学厚度（AOD）进行对比，从多源证据层面共同验证了本数据集在刻画咸海地区风蚀动态方面的可信度。本数据集可作为咸海地区土地退化评估、生态修复工程成效分析以及"绿色丝绸之路"生态安全研究的关键数据基础，同时也可为地理学、生态学与水土保持科学等领域的研究与教学提供数据支撑。
关键词：多源遥感；RWEQ模型；咸海；年际变化；中亚
DOI: https://doi.org/10.3974/geodp.2026.01.10
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2026.01.10
数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：
https://doi.org/10.3974/geodb.2025.12.01.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2025.12.01.V1.
咸海位于亚欧大陆腹地，曾是世界第四大湖泊，其生态系统的剧变是20世纪以来最严重的环境灾难之一。自1960年代以来，在气候变化与高强度人类活动的共同驱动下，咸海流域水资源平衡被彻底打破，湖体持续萎缩，至21世纪初已丧失超过90%的水域面积[1][2]。伴随湖床的大规模裸露，区域地表反照率、蒸散格局及局地气候条件发生显著改变[3][4]，进一步加剧了土地退化进程。干涸湖底沉积物在风力作用下成为盐尘暴与沙尘暴的重要源区，释放的大量含盐粉尘随西风环流长距离输送，对中国西北地区的整体生态环境构成潜在威胁[3]。咸海生态危机[5][6][7]不仅是中亚地区的区域性挑战，也是影响"绿色丝绸之路经济带"生态安全格局的重大环境问题，受到中国与中亚各国的高度关注。
在此背景下，咸海地区已成为研究干旱区土地利用-水文-气候耦合作用及其生态响应的典型区域[8]。由于土壤风蚀是引发咸海生态危机的关键驱动过程，系统开展该区域风蚀动态监测与风险评估研究，对落实中国—中亚峰会关于"推动解决咸海生态危机"的倡议、共建中国—中亚命运共同体具有重要现实意义。为支持相关研究与决策，高分辨率、长时序的土壤风蚀数据成为迫切需求。目前，遥感与风蚀模型的结合已被广泛应用于区域尺度风蚀评估[9][10][11]，其中修订版风蚀方程（Revised Wind Erosion Equation, RWEQ）因其参数可获取性强、物理机制明确，在干旱半干旱地区表现出良好适用性[12]。然而，现有数据产品在咸海周边，仍存在时空分辨率不足或验证薄弱的问题[13]。
本数据集基于多源遥感与再分析资料，在Google Earth Engine平台实现RWEQ模型的区域化应用与长时序模拟，可用于分析咸海萎缩过程中土壤风蚀的时空分异规律[14]，识别风蚀热点与演变趋势，并为区域土地退化评估、生态工程效益监测[15]及"一带一路"生态安全协同治理提供关键数据支撑。
《中亚咸海及周边区域风蚀模量数据集（1990–2020）》[16]名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 中亚咸海及周边区域风蚀模量数据集（1990–2020） |
| 数据集短名 | ASSR_WEMD_1990-2020 |
| 作者信息 | 余瑶，中国科学院新疆生态与地理研究所干旱区生态与资源科学数据中心，yuyao@ms.xjb.ac.cn 姚峰，中国科学院新疆生态与地理研究所干旱区生态安全与可持续发展国家重点实验室，yaofeng@ms.xjb.ac.cn |
| 地理区域 | 咸海盆地（北纬44°–47°、东经58°–62°） |
| 数据年代 | 1990–2020年 |
| 时间分辨率 | 年 |
| 空间分辨率 | 1 km |
| 数据格式 | .tif |
| 数据量 | 260 MB |
| 数据集组成 | 中亚咸海及周边区域风蚀模量数据 |
| 基金项目 | 新疆维吾尔自治区科学技术厅（2024E02030，PT2406） |
| 数据计算环境 | Google Earth Engine |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所 |
| 数据共享政策 | （1）"数据"以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用"数据"需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）"数据"的用户需要与《全球变化数据学报》（中英文）编辑部签署书面协议，获得许可；（4）摘取"数据"中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[17] |
| 数据和论文检索系统 | DOI，CSTR，Crossref，DCI，CSCD，CNKI，SciEngine，WDS，GEOSS，PubScholar，CKRSC |
本数据集集成了多源遥感、再分析与观测数据（表2），用于研究区域风蚀模量的计算与验证。
| 序号 | 数据类型 | 数据来源 | 数据内容 | 空间分辨率 | 时间范围 |
|---|---|---|---|---|---|
| 1 | 地面观测 | NOAA | 气象站数据（风速等） | – | 1990–2023年 |
| 2 | 气象数据 | ECMWF ERA5-Land | 风速、温度、降水、积雪深度 | 1 km | 1990–2023年 |
| 3 | 遥感 | MODIS | 地表温度、植被覆盖度、气溶胶光学厚度 | 1 km | 1990–2023年 |
| 4 | 土地覆盖数据 | ESA CCI | 全球土地覆盖数据 | 300 m | 1990–2023年 |
| 5 | 地形数据 | SRTM | 高程数据 | 30 m | – |
| 6 | 土壤数据 | HWSD | 土壤特性（质地、有机质等） | 1 km | – |
| 7 | 野外采样 | 哈萨克斯坦咸海干涸湖床 | 土壤采样数据（29个点位） | – | 2018年6月 |
| 8 | 粉尘/盐粒监测 | 咸海沙尘监测站 | 沙尘/盐分沉积数据 | – | 2000–2005年 |
本数据集采用欧洲中期天气预报中心（European Centre for Medium-Range Weather Forecasts，ECMWF）发布的ERA5-Land再分析数据集[1]作为主要气象数据源。该数据集提供了空间分辨率为0.1°（约10 km）的逐小时近地表风速、温度及降水估算值。为提升风速数据的可靠性，我们引入美国国家海洋和大气管理局（National Oceanic and Atmospheric Administration，NOAA）地面站涉及本研究区的日观测数据[2]，采用分位数映射法对ERA5-Land风速进行了系统性偏差校正。经验证，校正后的数据与实测值吻合度较高，决定系数（R2）超过0.87。
遥感数据集构成了空间分析的核心。本数据集采用中分辨率成像光谱仪（Moderate-resolution Imaging Spectroradiometer，MODIS）数据提取1-km空间分辨率的归一化植被指数（Normalized Difference Vegetation Index，NDVI）、气溶胶光学厚度（Aerosol Optical Depth，AOD）与地表温度（Land Surface Temperature，LST）。植被指数来源于美国地质调查局/美国国家航空航天局MOD13A2产品[3]，并通过Z-score标准化处理以降低传感器特异性偏差。土地覆盖数据获取自欧洲航天局气候变化倡议项目300-m分辨率数据集[4]，采用众数重采样方法聚合至1-km分辨率。地形数据[5]通过30-m分辨率的航天飞机雷达地形测绘任务获取，采用均值重采样方法升尺度至1 km以匹配其他数据源。核心的MODIS数据均经过严格质量控制，包括基于QA波段（置信度阈值>90%）的云掩膜处理，并对NDVI等参数开展了时序一致性检验（剔除偏离五年移动平均值±3σ的异常值）。
为保障时序一致性与年际可比性，所有数据集均统一重投影至WGS 1984地理坐标系，并采用双线性插值法（连续变量）与最近邻法（分类数据）重采样至统一的1-km空间分辨率。针对NDVI、AOD及降水数据，采用五年移动平均法进行时序平滑处理，以有效抑制短期异常波动与传感器噪声的影响。
土壤属性信息提取自空间分辨率为1 km的世界和谐土壤数据库[6]，包含土壤质地、有机质及碳酸钙含量等关键参数，这些是计算RWEQ模型中土壤可蚀性因子的核心输入项。通过对29个野外采样点数据进行验证，土壤参数误差严格控制在容许范围内（质地±5%，有机质±0.5%）。
修正风蚀方程（RWEQ）是一个基于经验发展而来的风蚀估算模型。其核心机理在于综合模拟气候（如风速、降水）、土壤（如可蚀性成分含量）及地表覆盖（如植被、糙度）等多要素对风蚀过程的共同影响，现已成为量化干旱与半干旱地区土壤风蚀强度的常用工具[18][19][20]。
在本研究中，RWEQ模型的关键输入参数来自遥感和气象数据集；单位面积的风蚀速率计算如下：
(1)
(2)
(3)
式中，
为土壤侵蚀速率（t/(hm2·a)），
为最大土壤传输速率（kg/m），z为最大侵蚀发生的下风距离（m）。模型假设最大风蚀发生在田地中点（z ≈ S/2），S为地块长（m），
为地表粗糙度因子，WF为气候因子，EF为土壤可蚀性因子，SCF为土壤结皮因子，C为植被覆盖因子。
气候因子（WF）的计算公式为:
(4)
式中，SW表示土壤湿度（%），SD为积雪覆盖因子，u2代表2-m高度处实测风速（m/s），u1是2-m高度的阈值风速（假设为5 m/s）（u2>u1），N为风速观测时段数，Nd指天数（d），
为空气密度（kg/m3），g表示重力加速度（m/s2）。
土壤可蚀性因子（EF）与土壤结皮因子（SCF）由以下公式给出：
(5)
(6)
式中，SA、CL和SI分别代表砂土、黏土和粉土含量（%），OM表示土壤有机质含量（%），而CaCO3代表碳酸钙含量（%）。植被覆盖因子（C）针对5类土地覆盖类型（森林、灌木地、草地、耕地和裸地）进行计算，公式如下：
(7)
式中，
为植被特定系数，SC表示基于NDVI得出的植被覆盖度（%）。这些植被特定系数引自Fryrear[21]针对干旱地区的研究，取值分别为0.153,5（森林）、0.092,1（灌木地）、0.151,1（草地）、0.043,8（耕地）和0.076,8（裸地），与研究区土地覆盖类型相符。
本数据集采用谷歌地球引擎（GEE）作为计算平台，其独特的云端架构有效解决了大范围风蚀模拟中的3大关键问题：首先，平台内置的并行计算能力支持对咸海地区1990–2020年共31年时序数据的高效处理；其次，GEE直接集成了多源地理空间数据（包括MODIS地表反射率产品、ERA5-Land气象再分析数据和SRTM高程数据），实现了数据的无缝调用与统一管理；第三，其优化的空间分析功能在保持1-km分辨率像元级计算精度的同时，确保了区域尺度模拟的可行性。基于GEE的集成化框架展现出卓越的计算效率，完整生成31年风蚀模量年际空间数据集仅需1天时间，为长时序、大范围风蚀模拟研究提供了关键技术支持。
本研究的技术流程如图1所示，该计算框架包含3个主要步骤：首先，在数据预处理阶段，气象数据（ERA5、NOAA）和土地覆盖数据（MODIS、ESACCI）被重采样至1 km分辨率以确保空间一致性，并通过偏差校正技术修正风速数据以提高精度；其次，应用RWEQ模型方程计算1990至2020年间风蚀速率的时空分布；最后，利用优化后的模型模拟1990–2020年咸海流域月尺度土壤风蚀潜力，生成连续时空序列数据集，并利用模拟的侵蚀速率与沙尘通量观测数据及卫星反演的气溶胶光学厚度（AOD）数据进行验证和校准。
图1 数据研发技术流程图
下载中亚咸海及周边区域风蚀模量数据集包含1990–2020年共31期土壤风蚀模量空间数据，以.tif格式存储，空间分辨率为1 km，单位为kg/m2/y。数据文件命名规则为：ASSR-WEMD-年份.tif，其中ASSR代表"中亚咸海及周边区域"，WEMD为"土壤风蚀模量"，年份为四位数纪年。例如，ASSR-WEMD-2010.tif表示2010年该区域的土壤风蚀模量栅格数据。
本数据集所有图层均采用统一的地理坐标系与像元大小，便于进行时序分析与空间建模。每个.tif文件包含1个单波段浮点型栅格层，像元值即表示该位置的年均风蚀模量。为体现数据集的典型样本特征，选取以下3个代表性年份数据进行说明：（1）ASSR-WEMD-2005.tif，代表序列中风蚀模量最低的年份，反映早期至中期风蚀活动的减弱阶段；（2）ASSR-WEMD-2010.tif，反映风蚀模量发生显著跃升的关键转折年份；（3）ASSR-WEMD-2015.tif，对应整个时序中风蚀模量达到峰值的年份，标志侵蚀强度最高阶段。
图2展示了1990–2020年共31个年度的风蚀模量空间分布。从整体年际变化来看，咸海盆地尤其是东部干涸湖床区域始终是侵蚀的核心区，且侵蚀强度总体呈加剧趋势。然而，自2012年起出现了一个重要转变，即湖床内部的侵蚀有所减弱，而周边区域的侵蚀强度明显上升，这表明侵蚀范围正在扩大，且主导侵蚀源已从裸露的湖床中心逐渐迁移至植被退化、土壤脆化的湖岸带及毗邻区域。
图2 土壤风蚀模数年际空间分布图（1990–2020）
下载在1990–2020年整个研究期间，全区土壤风蚀模数以年均0.33 kg/m2的速率显著递增，空间变异率介于–1.2至5.09 kg/m2。其中，侵蚀剧烈加剧的区域（年增长率>1.8 kg/m2）主要分布在以下3类地区：咸海中部新近干涸的湖床、克孜勒库姆沙漠的活跃沙区，以及西部高原的脆弱地表区。这些区域共同特征是地表物质松散、植被覆盖极低，因而对风力作用和环境变化响应最为敏感。与之形成鲜明对比的是，咸海北部残留水体周边及锡尔河、阿姆河下游三角洲的灌溉绿洲区，年变化率极低（<0.1 kg/m2），这主要得益于水体或灌溉维持的较稳定植被覆盖，有效抑制了风蚀发展。
为评估1990–2020年中亚咸海地区风蚀模量年际空间数据集的准确性与可靠性，本研究将模拟结果与咸海沙尘监测站（2000–2005年）记录的粉尘通量实测数据进行对比验证，发现模拟风蚀强度空间格局与实地观测数据具有较高一致性。
首先，利用咸海地区监测站（2000–2005年）记录的实测粉尘通量数据与对应位置的模拟风蚀模数进行对比分析，结果显示两者在空间格局上高度吻合，且具有强相关性（r=0.72, p<0.05），表明模型能够可靠捕捉区域风蚀的时空变化特征。进一步，我们将模拟的年际风蚀模量与卫星反演的气溶胶光学厚度（AOD）数据进行区域关联分析，发现二者在年际尺度上亦呈高度显著相关（r=0.85, p<0.001），这说明模拟结果不仅与地面观测一致，也与反映大气粉尘负载的遥感指标相吻合，从而从多源证据层面共同验证了本数据集在刻画咸海地区风蚀动态方面的准确性与可靠性。
本数据集以.tif格式存储，用户可在ArcGIS、QGIS、ENVI、Google Earth Engine等主流地理信息与遥感平台中直接读取，并进行可视化、空间查询、统计分析与制图输出。数据适用于咸海区域1990–2020年间土壤风蚀模数的年际变化趋势分析、空间格局演变识别及驱动机制研究，也可作为风蚀风险评估、土地退化监测和生态治理效果评价的基础数据。此外，本数据集支持作为模型输入数据用于未来风蚀情景模拟，为咸海生态修复与区域可持续发展决策提供科学依据。
本研究及未来工作存在以下可改进与深化之处：（1）模型与输入数据的优化，当前模型主要依赖卫星遥感与再分析数据驱动。未来可融合更高时空分辨率的遥感产品（如Sentinel系列数据）、增加实地观测的风速及土壤属性等关键参数，以进一步提升模型的模拟精度，尤其是在表征局部微地形与土壤异质性方面。（2）过程机理的精细化表征，现有框架主要基于修订版风蚀方程（RWEQ）进行年际尺度的模拟。未来可集成更复杂的风蚀过程模型，或耦合水文、植被动态模型，以更精细地刻画日/季节尺度风蚀事件、土壤湿度与冻融过程对风蚀的影响，以及植被动态与风蚀的反馈机制。（3）应用场景的拓展与决策支持，未来可将本数据集与区域气候模型、土地利用变化情景相结合，开展高分辨率的风蚀风险预测预警研究。同时，可探索将风蚀模量数据与生态系统服务评估模型、社会经济数据相集成，量化风蚀灾害的经济社会影响，从而为制定更具针对性的生态修复与适应性管理策略提供更强大的决策支持。
通过对上述方面的持续改进，有望进一步深化对咸海地区风蚀过程的理解，并为全球其他受风蚀威胁的干旱半干旱区提供更精准、更具前瞻性的评估工具与数据基础。
作者分工：余瑶承担数据统计工作，并撰写了数据论文；姚峰对数据集的开发做了总体设计、数据质量控制和论文修改；李昌俊承担数据验证工作。
利益冲突声明：本研究不存在研究者以及与公开研究成果有关的利益冲突。
[1] ERA5-Land. https://cds.climate.copernicus.eu/datasets/reanalysis-era5-land?tab=documentation.
[2] NOAA日观测数据. https://www.ncei.noaa.gov/data/global-summary-of-the-day/archive/.
[3] MOD13A2. https://ladsweb.modaps.eosdis.nasa.gov/.
[4] ESA CCI. https://maps.elie.ucl.ac.be/CCI/viewer/.
[5] SRTM. https://dwtkns.com/srtm30m/.
[6] HWSD. https://gaez.fao.org/pages/hwsd.